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This work explores the potential reach of the 7 TeV LHC to new colored states 
in the context of simplified models and addresses the issue of which search regions 
are necessary to cover an extensive set of event topologies and kinematic regimes. 
This article demonstrates that if searches are designed to focus on specific regions of 
phase space, then new physics may be missed if it lies in unexpected corners. Simple 
multiregion search strategies can be designed to cover all of kinematic possibilities. A 
set of benchmark models are created that cover the qualitatively different signatures 
and a benchmark multiregion search strategy is presented that covers these models. 
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There is a place where the sidewalk ends 
And before the street begins, 
And there the grass grows soft and white, 
And there the sun burns crimson bright, 
And there the moon-bird rests from his flight 
To cool in the peppermint wind. 

- "Where The Sidewalk Ends" Shel Silverstein, 1974 

1. INTRODUCTION 

By the end of 2011, the LHC is expected to have accumulated C(lfb _1 ) of integrated 
luminosity at \fs = 7TeV. At this early stage, models of new physics with large production 
cross sections should be prioritized in experimental analyses. Among those are spectra with 
colored particles that decay into jets and a stable invisible particle, which are a common 
prediction of extensions of the Standard Model motivated by the hierarchy problem and the 
dark matter puzzle. 

The Tevatron carried out a variety of searches in jets and missing energy and extended 
bounds on colored objects in specific contexts [IrS]. Most searches, however, were optimized 
for mSUGRA-type benchmarked scenarios, that are affected by strong assumptions on the 
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spectrum, mass splittings and branching ratios, and therefore underrepresent the kinematic 
possibilities and decay topologies. At the LHC, both on the theory and experiment sides, 
similar model-specific studies for new physics prospects have been carried out [5TI2T] . More- 
over, many analyses are often obscured by the presentation of the results in terms of high 
energy mSUGRA parameters such as m and m^, making it non-trivial to translate the 
bounds for alternative theories. Previous work recast the Tevatron bounds for more general 
scenarios and showed that its reach could have been significantly extended, had a reanaly- 
sis been performed with a less benchmark driven and more comprehensive search strategy 
[22j [23] . Some model- independent searches have been carried at the Tevatron [21H26] . but 
the inferred limits from these searches have never been performed and are difficult to do a 
posteriori. 

An alternative paradigm for creating searches and exploring common features of new 
physics was recently put forward and dubbed "simplified models." Simplified models pa- 
rameterize the new physics by a simple particle spectrum, its production mode and decay 
topologies with the masses, cross sections and branching ratios taken as free parameters 
[22] |23j I2TH3T] . These simplified models capture generic kinematic properties of models that 
are relevant for early searches. Particles that are not involved in a specific signature are 
decoupled from the simplified model. In some sense simplified models define a perturba- 
tion series in model space that approximates specific theories and facilitates their mapping 
into experimental observables. While still model-dependent, simplified models help reduce 
model dependence that often plague top-down parameterizations of new physics. Simplified 
models have been used in the context of solving the "LHC Inverse Problem" [22] and also in 
characterizing experimental anomalies [28J. Other approaches in reducing model dependence 
of searches have been to more fully explore the 19 parameter "phenomenological MSSM" 

una Esi ■ 

In this article, the early LHC discovery potential for jets and missing energy signatures 
is studied in the framework of simplified models. The simplified spectrum considered here 
consists of a gluino-like object that is pair-produced and decays to jets and a stable neutral 
particle that escapes detection. The decay can proceed directly to the neutral particle or 
through a cascade, in which case the spectrum is augmented by intermediate particles. The 
decay modes were chosen so as to cover a diverse kinematic range and correspond to two- 
body direct decays, three-body direct decays, three types of one-step cascades and a two-step 
cascade. The scope of these simplified models is broader than the realm of supersymmetric 
theories [351410] as long as spin correlations are irrelevant in the discovery process [JT]. 

The first result of this article is a systematic quantification of the reach for these colored 
objects with an extensive set of combined cuts on (missing energy) and (visible 
energy). These results are used to subsequently extract the optimal sensitivity of early LHC 
searches to the gluino branching ratios times production cross section. For a QCD reference 
cross section and 100% gluino branching ratio to jets plus $j>, the LHC may find a 2a 
evidence for gluinos as heavy as 800 GeV in models with large mass splittings, with lfb -1 
of integrated luminosity. For compressed spectra, the estimated current limit of 150 GeV 
on gluino masses can be extended up to 500 GeV. 

The optimal reach for the set of simplified models considered is not obtained by a single 
set of cuts on any set of variables, but instead requires searching in several distinct regions. 
The second result of this article presents a minimal set of search regions designed to obtain 
a near-optimal reach for each simplified model. Six search regions are necessary to cover the 
entire mass parameter space of the simplified models, guaranteeing that signatures of new 
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physics in all hadronic jets + MET channels will not be missed. These six searches are: 

• dijet high MET for compressed spectra 

• trijet high MET for heavy gluinos decaying via 2-body direct decay to the LSP 

• multijet low MET for light gluinos decaying via 3-body direct decay or cascade 
decay to the LSP 

• multijet moderate MET for intermediate mass gluinos decaying via 3-body direct 
decay or cascade decay to the LSP 

• multijet high MET for heavy gluinos decaying via 3-body direct decay or cascade 
decay to the LSP 

• multijet high H T for heavy gluinos decaying to light LSP's. 

The third result of this article is a short list of benchmark simplified models that can be 
used to ensure that all regions of parameter space are being covered with equal diligence. 
These benchmark models are useful because they provide examples of spectra that are not 
represented in mSUGRA or other existing parameterizations. 

The outline of this article is as follows. Sec. [2]provides a precise definition of the simplified 
models used in this study. Sec. [3] describes the Monte Carlo (MC) implementation of signal 
and background. These MC calculations are used on Sec. [4] for the estimation of the cross 
section sensitivities, and on Sec. [5] for investigating and designing a comprehensive search 
strategy on $t and Ht- Sec. [^explores how to use the results of single topology simplified 
models to multiple topology simplified models. The conclusions are summarized on Sec. [7l 
Finally, App. |a| shows the estimated reach for the simplified models with 45 pb _1 and 1 fb~ 
and App. [B] gives a set of fully specified benchmark simplified models. 

2. SIMPLIFIED MODELS FOR COLORED OCTETS 

The simplified models described in this section are effective field theories after electroweak 
symmetry breaking. Therefore the states of the theory will have well-defined quantum 
numbers under SU(3) C x £7(1) EM . The two principle states in the set of simplified models 
are a color octet Majorana fermion, g, and a stable neutral Majorana fermion, x°- The 
topologies considered in this article will all commence with g being pair produced through 
its QCD interactions. The fully differential production cross section of g is assumed to be 
proportional to the fully differential tree-level QCD production cross section. This is a good 
approximation whenever non-Standard Model states do not play a dominant role in the 
production of g. Important cases such as resonant production are not directly proportional 
to the differential QCD production cross section and deserve a separate treatment. The 
results derived in Sec. [4] and Sec. [5] are parametrized in terms of m- g and m x o and the 
product of the inclusive production cross section of g and any relevant branching ratios. 
By allowing the total g production to be rescaled, many different theories are effectively 
described by a single simplified model. Effects such as different g multiplicities or spin 
degrees-of-freedom and even t-channel squarks in supersymmetric theories are well-captured 
by this parameterization of simplified models. For convenience, the inclusive production 
rates will be referenced against the QCD next-to-leading order cross section. 
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The decay of each g into x° an d additional jets proceeds through an effective operator. 
The most straight-forward augmentation of this simplified model is to add intermediate 
particles in the decay chain of g which can significantly alter the decay kinematics. Two 
benchmark particle types are chosen in this study: a charged Dirac fermion x ± with the 
quantum numbers of a chargino and a neutral Majorana fermion x'° with the quantum 
numbers of a neutralino. In addition to the direct decay of 

9 -> X° + X 

where X = g or qq, two other decay routes are studied: a one-step cascade decay, 

~g -> X ± + X -+ X ° + W ± +X, 
and a two step cascade decay, 

g ->■ x ± +X ->■ X '° + W ± + X -> x° + ^° + + X 

In the cascade decays of the g down to x°, leptons may appear as by-products of the 
W ± and Z° decays. Although leptons may be present in the final state, this study will be 
concerned exclusively with all hadronic final states plus missing energy. Leptonic channels 
can in principle increase the sensitivity to simplified models with electroweak gauge bosons in 
the final state; however, these decay topologies are "lepton poor" due to the relatively small 
branching fraction of electroweak vector bosons into leptons. In this study, the sensitivity 
to g one-step decays through a chargino in the one-lepton channel was found to be worse 
than that in the all hadronic mode. 

The four g decay modes considered are detailed below and illustrated in Fig. [11 

2.1. Two- body direct decay 

9 9X° 

In this simplified model, the dominant decay mode is g — > x° + 9- This decay proceeds 
through an effective operator of the form 

Ant = ^~ga^ X °G^ + h.c. (1) 

This decay mode can occur in gauge mediated supersymmetry breaking models when the 
gluino is the next-to-LSP and %° is the gravitino. Two body g decays to x° an d a gluon also 
appear at loop level in standard supersymmetric theories and can become large when there 
are sizeable splittings between the left and right handed squarks and can even become the 
dominant decay mode in split supersymmetric models [221 123] • 

In Universal Extra-Dimensions (UED), the KK-excitation of the gluon may dominantly 
decay to a gluon and a KK-graviton [35]. Note that in this case the pair-produced octet is a 
massive spin-one particle and the invisible final state is a massive spin-2 particle. Neverthe- 
less, such UED scenarios can still be parametrized by the simplified model considered here 
given that the inclusive production cross section is taken to be a free parameter and most 
of the signature's features are determined by kinematics alone. 

In addition to theories where the colored particle is a color octet, this simplified model also 
well-approximates theories where the new colored particle lies in a different representation 
but have the same decay topology, such as squark NLSPs. The approximation breaks down 
when initial or final state radiation becomes important. 
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FIG. 1: Simplified spectra and decay modes considered in this study. The masses of the pair 
produced colored state, g, and of the the invisible final state, x , are varied over the range from 
100 GeV to 1 TeV. The masses of the intermediate states in the decay chains are held at a fixed 
distance from m,g and m^o, corresponding to slices through the mass parameter space that are 
representative of all the kinematical configurations. 



2.2. Three-body direct decay 

g ->• qq'x 

In contrast to the previous simplified model, the gluino decay may still be direct, but 
proceed via 3-body. In supersymmetric models with decoupled squarks, e.g. in split- 
supersymmetry [44J, the gluino can decay to an electroweak gaugino and two light flavored 
quarks. This decay proceeds through a dimension-6 operator 

Ant = -^kgqiQiX + h.c. . (2) 

where i runs over the different quark flavors and cji X is the Yukawa coupling between the 
quark-squark and and x°- This article is assuming that the decays proceed into light fla- 
vored quarks. The final state flavor structure is determined by the mass spectrum of the 
corresponding squarks with the decay through lighter mass squarks occurring more rapidly. 
Assuming all electroweakinos are kinematically accessible in gluino decays, direct three body 
decays dominate in the following cases 

• %° = B and the right handed squarks are lightest, 
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• X° — W and the left handed squarks are lightest, 

• x° — H and the heavy flavor quarks are kinematically accessible in gluino decays . 

In AMSB scenarios, the LSP is usually wino-like and there no strong hierarchy in squark 
masses; however, due to the fact that the wino gauge- Yukawa coupling is larger, the direct 
decay of the gluino has a large branching ratio. In mSUGRA and GMSB-like models, the 
LSP is usually bino-like and there is no strong splitting between the left and right-handed 
squarks; therefore, the direct decays usually do not dominate. The study of heavy flavor 
decays are beyond the scope of this article, see [35] . 

This simplified model is relevant not only for 3-body direct decays, but it also effectively 
describes cascade decays in which there is a large mass splitting between the gluino and the 
first intermediate state in the cascade and the rest of intermediate states are compressed 
near the LSP. This occurs frequently when the LSP is an SU(2)l multiplet such as the wino 
or Higgsino. In these scenarios, the mass splittings between the LSP and NLSP are typically 
less than 50GeV and can be as small as (9(100 MeV). 

2.3. One-step cascade decay 

g ->• qq'x ± -> qq'{W ± x°) 

In this simplified model, the gluino goes through a 3-body direct decay to a chargino 
that subsequently decays to a gauge boson and the LSP. This simplified model is commonly 
realized in mSUGRA, and more generally the chain gluino — > heavy electroweakino — > 
lightest electroweakino is preferred in many supersymmetric scenarios [36] . A similar chain 
KK-gluon — > KK-gauge boson — > KK-graviton is also present in Extra Dimensions, because 
spin correlation are a mild effect on early discovery. 

One of the challenges in non-minimal simplified models is the proliferation of parameters. 
There are two additional parameters to consider in one-step cascade decays: the mass of the 
intermediate particle m x ± and the branching ratio of g decaying into x ± - Each simplified 
model will initially be considered with branching ratios set to 100%. Multiple decay modes 
can be studied by taking linear combinations of single decay modes. This is discussed in 
Sec. El 

By setting branching ratios to 100%, the number of parameters in one step cascades are 
reduced to four, namely m§, m x o, m x ±, and a(pp —> gg + X). The choice of m x ± alters 
the kinematics of the theory and therefore all four parameters are important. This article 
explores this four dimensional parameter space by choosing several distinct relationships 
between the three masses of the theory (or "mass slices"). The mass slices were chosen to 
maximize the distinctive features of one-step cascades and capture all the relevant corners 
of phase space. Those are specified by the following choice of intermediate chargino masses: 

m x ± = m x o + r{rrig — m x a), (3) 

where three values of the parameter r chosen are: r =1/4, 1/2 and 3/4. The case of r = is 
identical to the direct 3-body decay. The case of r = 1 closely resembles the direct two-body 
decay when the boost of the W ± becomes large and its decay products merge together. 
For the mass ranges of interest for the LHC, the W ± is a light particle and therefore most 
theories are in this boosted W ± regime. 
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An alternative mass slice that could be studied is 

m x ± ~ m x o + m w ± (4) 

in order to explore the effects of on-shell decays near threshold. This mass slicing is present 
in a subspace of the three parameterizations adopted in Eq. [3] when 

m x o = m g — m w ± jr. (5) 

Threshold effects are fairly modest because the mass scales accessible at the LHC are suf- 
ficiently above mw±. The results presented in Sec. [4] and App. [A] do not show anomalous 
behavior near this line. Threshold effects are important for lighter g masses. In |31] , this 
can be seen as a sharp drop in the cross section sensitivity along the line in Eq. |4} 

Given that the intermediate particle in this simplified model is a chargino, all events 
have two W ± bosons in the final state. Clearly, alternative simplified models resembling 
this one exist, in which the intermediate state is neutral and decays to a Z° boson instead of 
a W ± . One of the open questions in the study of simplified models is how adjacent theories 
in "model space" approximate each other. This example of exchanging a neutralino for a 
chargino, or equivalently a final state W for a Z°, is a prime case study for gaining intuition 
for the scope of simplified models. 

When exchanging a W ± for a Z°, the mass difference is a small effect at the LHC. 
The primary difference is the presence of leptonic events which get vetoed in the searches 
considered in this study. Ultimately, there are two questions that need to be addressed. 
The first is whether the optimization of search regions is affected by the choice of final 
state vector bosons. This question is the most critical because it can lead to a delay in the 
discovery of new physics. Fortunately, the search design is insensitive to the small difference 
in the number of leptonic events between the two theories. 

The second question is how to translate the limits from one simplified model to the other. 
Answering this question requires understanding the differences in the acceptances/efficiencies 
for events with Z°-final states versus W /± -final states. The efficiencies are similar for most 
mass spectra, with larger discrepancies arising in regions of heavy g and light x°-> differing 
by typically 20% due to the 20 GeV lepton veto. This results in a slight gain in sensitivity 
for simplified models with Z° bosons in the final-state. To understand the magnitude of 
this difference, consider the following back-of-the-envelope calculation. The lepton veto is 
not perfectly efficient in identifying leptons. The efficiency for finding isolated individual 
leptons in the jet-rich environment of g decays is estimated by PGS4 to be 

e e ~ 74% (6) 

with approximately equal efficiencies for both electrons and muons. This causes a leakage of 
leptonic events into the signal region. The events in this sample have two W bosons with 
branching ratios 

B 2Wm = Br (2W -> hadrons) ~ B 2 Wh ~ 62% 

B 2Wu = Br (2W ->■ £ + hadrons) ~ 2B Wh (l - B Wh ) ~ 33.5% 

B 2W2e = Br (2W -> 21) ~ (1 - B Wh f ~ 4.5%. 

where 



B Wh = Br (W qq) + Br (W rv) ~ 82.4%. 
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This results in a leakage fraction of leptonic events into the all-hadronic channel of 

L 2W = (1 - e t )B 2Wu + (1 - e e ) 2 B 2W2e ~ 9%. (7) 

Thus the resulting fraction of 2W ± events passing the lepton veto is 

e 2 ^ eto ~ B 2Woe + L 2W ~ 71%. (8) 

The analogous reasoning for events with 2Z° in the final state leads to 

e 2 e z veto ~ 88% e£ veto /e 2 ^ eto - !- 2 3. (9) 

This back-of-the-envelope calculation is surprisingly accurate for the large mass splitting 
regions and accounts for the 20% difference in sensitivity mentioned above for W versus 
Z° simplified models. For more compressed spectra where the W ± /Z 0l s do not have phase 
space to go on-shell, the final-state leptons are so soft that the leptonic "leakage" rates into 
the hadronic region can easily go up to 100%. In such cases there is virtually no distinction 
between the two types of simplified models, implying nearly identical search sensitivities. 

2.4. Two-step cascade decay 

9 W'X ± "> qfiW+x' ) -> qq / (W ± (Z\ )) 

The last simplified model in this study consists of g decaying through two intermediate 
states to v°. These are encountered in the MSSM if the ordering of the electroweakinos is 
W — H — B or B — H — W because the decay widths of W <-> B are suppressed relative to the 
B <r¥ H -H- W transitions. In extensions to the MSSM with singlets such as the NMSSM, 
the decay chain g — > W — > B — > S is a common cascade because the singlino, S, can have a 
small coupling to the MSSM. 

Higher-step cascades are possible, but they typically have small branching ratios due to 
phase space suppression and the sheer number of combinatoric possibilities. Furthermore, 
the kinematics frequently resemble the simplified models already considered. 

The number of parameters in this two-step cascade grows considerably over the three in 
the direct decay simplified models and five in one step cascade decays. In two-step cascade 
decays there are four masses, three branching ratios and the production cross section. As 
with the one-step cascade decay, the branching ratio for the two-step cascade decay is set 
to 100%. Sec. [6] considers going away from this limit. 

In order to simplify the four dimensional mass parameter space, the following mass slicing 
is chosen for the intermediate states 

m x ± = m x o + r(rn~ g — m x o) 

m x io = m x o + r'(m x ± — m x o), (10) 

with r = r' = 1/2. More choices could be studied, but this specific parameterization 
maximizes the difference in the kinematics between the two-step cascade decay and the 
one-step and direct decay simplified models previously considered. 

The effects of the W ± and Z° going on-shell in the cascade decay occur for the following 
values 

m x o = nig — m z o/rr' Z° on-shell threshold 

m x o = nig — m w ±/r(l — r') W ± on-shell threshold. (11) 
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Only one of these can be satisfied at a time unless r' = 1/(1 + mw±/mzo) — 0.53 which is 
close to the value chosen. The threshold effect can be seen in the results in Sec. H] and in 
App. [A] as a loss of sensitivity near the line m- g ~ m x o + 350 GeV. 

This simplified model has g decaying into a charged x ± - subsequently decays to a 
W boson and a neutral x'°- Finally, x'° decays to a Z° boson and x°- The final state 
therefore contains four jets, two W's, two Z 0> s, and two x°- Alternate charge choices for 
the intermediate states can be made and result in different combinations of vector bosons 
in the final state such as those containing AW ± or AZ°. As in the one-step cascade decay, 
the lepton veto is responsible for the dominant differences between the sensitivity to these 
final states. A similar estimation can be made for the differences between 2W ± 2Z° final 
states versus AW ± or 4Z° final states. The resulting search sensitivity to 2W ± 2Z° simplified 
models is typically 16% worse than for 4Z° simplified models and 16% better than for AW ± 
simplified models. 



3. BACKGROUNDS AND SIGNAL SIMULATION 

The dominant Standard Model backgrounds to jets and $t signatures are tt + jets, 
W^-l- jets, Z°+ jets, and QCD. In this study, subdominant processes such as single top and 
diboson production in association with jets are not included. The matrix elements for parton 
level events are computed in MadGraph 4.4.32 [U] with CTEQ6L1 parton distribution 
functions used throughout |48j. Variable renormalization and factorization scales are set to 
the transverse energy of the event [49J. The SM parton level processes generated are 

pp -> W ± + nj l<n<3 (12) 

pp Z° + nj 1 < n < 3 

pp tt + nj < n < 2 

pp -> W ± + bb + nj < n < 2 

pp-> Z° + bb + nj < n < 2 

where j stands for gluons and light flavored quarks. The Z° and W plus jets samples are 
forced to decay to final states involving neutrinos, but no restrictions on tt events are placed. 

The Standard Model contribution to missing energy distributions peaks at low energies, 
whereas many signatures of new physics yield events with extremely high missing energy. 
Therefore it is important to have sufficient Monte Carlo statistics on the tail of the $t 
distribution. To achieve sufficient statistics, different samples are generated for each SM 
process, where each sample has the pr of the massive particle lying in a given interval. For 
instance, Z° + jets parton level events are broken into three samples with 

GeV < p T;Z o < 200 GeV 
200 GeV < p T>z o < 300 GeV 
300 GeV < p TtZ o . 

In the case of tt, where there are two heavy particles, the samples are divided by the larger 
Pt of either top quark in the event. 

Contributions from QCD to jets and $t can come from either detector effects and jet 
energy mismeasurement, or neutrinos appearing in the decay of heavy flavor hadrons. To 
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FIG. 2: Schematic representation of the different samples generated for QCD 4j processes. The 
events in a given sample have the pr of the two hardest jets within a chosen range. The lower right 
box can not be populated the four jet events because Pt ^ > Pt j 2 + Pt j 3 + Pt j 4 - 



estimate the QCD contribution to jets and $t signatures the following subprocesses are 
generated 

PP ~> JJJJ 

pp —J- bb + nj < n < 2 (13) 

The QCD 4j background is calculated by subdividing the samples by the pt of the two 
leading jets in the event at parton level. Eight separate samples are generated, specified in 
Fig.|2} ' 

The signal, g pair-production, was generated in association with up to two jets at parton 
level, 

pp — »• gg + nj < n < 2. (14) 

The effects of including additional radiation in signal processes have been documented in 
several studies [22| [231150]. 

For both signal and backgrounds, the showering, hadronization, particle decays, and 
matching of parton showers to matrix element partons were done in PYTHIA 6 . 4 [21] . Match- 
ing of parton showers to matrix elements is an important step to correctly utilize the multiple 
final state samples generated above. The matrix elements better describe hard radiation, 
while the parton shower generates softer radiation that fills out jets [52J. The MLM parton 
shower/matrix element matching scheme is used with a shower-fc^ scheme introduced in [53J. 
The matching scales used in this article are listed in the table below. 



Sample 


Q Match 


tt + jets 


100 GeV 


V+jets 


40 GeV 


V + bb + jets 


40 GeV 


bb + jets 


50 GeV 


gg + jets 


100 GeV 



(15) 



Hard jets beyond the multiplicities listed in (12) must be generated by the parton shower. 
In particular, for W ± + jets and Z° + jets, the fourth jet and beyond are generated through 
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the parton shower. This approximation has been validated by several studies. For instance, 
in W + + jets, the discrepancy in the inclusive rate for four jets and $t from matching up 
to three jets versus four jets is 0(15%) [52"] . 

Next-to-leading-order (NLO) corrections alter the predictions of both signal and back- 
ground. With parton shower /matrix element matching, the shapes of differential distribu- 
tions are accurately described by tree level predictions. The largest corrections are to the 
inclusive production cross section and can be absorbed in i^-factors. The leading order cross 
sections of the signal are normalized to the NLO cross sections calculated in Prospino 2 . 
|54j . The ti + X leading order production cross section is scaled to the NLO one with the 
same K factor used in [29] obtained from (55]. As a cross-check, the tt+jets $t distribution 
is checked to agree with that from [56]. Similarly, the W/Z+ jets leading order cross sections 
are scaled with the same normalization from J5BHE0], and checked to agree with |5"6"] . 

PGS 4 is used as a "transfer function" that takes hadron level events to reconstructed, 
detector-level objects that represent the objects that experiments at the LHC report results 
on [BT] . This study uses the PGS 4 ATLAS card and this has been shown to reproduce 
results to 0(20%) accuracy. 

One of the drawbacks of PGS 4 is that it uses a cone jet algorithm with AR = 0.7. This 
is an infrared unsafe jet algorithm, but better represents the anti-A;^ algorithms used by the 
experiments than the algorithm. The Standard Model backgrounds change by at most 
(9(10%) when varying the cone size to AR = 0.4 . 

The signal offers a more varied testing ground for the effects of changing the jet algorithm. 
Two competing effects are found. The first is that there is more out-of-cone energy for smaller 
cones, resulting in less energetic jets. The second effect is that smaller cone jet algorithms 
find more jets. These two effects are more pronounced when contrasting a jet-poor 2-body 
decaying g with a jet-rich 2-step cascade decaying g. The dependence of the kinematic cut 
efficiencies on AR varies with mass splitting between the g and x°- F° r compressed spectra, 
when the pt of the jets is reduced, the efficiencies for the smaller cone size decrease because 
jets fall below the minimum jet pt requirement. For widely spaced spectra, where jets are 
energetic, more jets are found with a smaller cone size and the efficiency to have multiple 
jets passing the minimum jet requirement increases. Altogether, the efficiencies differ by at 
most 0(20%) and is consistent with other studies [H2] • This effect is not included in the 
remainder of this article, where the analysis is performed with a fixed AR = 0.7. 

4. OPTIMAL SENSITIVITIES 

One of the goals of this article is to explore the necessary kinematic cuts that will opti- 
mally distinguish signal from background. The simplified models described in Sec. [2] provide 
a varied signal space when the full range of rrig and m x o is considered. The first step is 
setting which kinematic variables will be investigated. Several options were explored: 

• missing transverse energy, $t, 

• visible transverse energy, Ht-, 

• transverse momentum of the leading jets, pTi with i = 1,2, 3,4, 
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• effective mass, M^, 1 

• fractional missing energy, $ T /M eS , 

• fractional jet momenta, e.g. ptiJ Pt ■ 

Combinations of these kinematic variables were considered in defining search regions in 
conjunction with varying multiplicities of jets. Since the primary concern of this article 
is the discovery of new states, as opposed to measurement of their properties, the most 
effective strategy is to keep the largest fraction of the signal possible. Placing cuts on several 
variables, as in a multivariate analysis, has the drawback that it increases the complexity 
of the search. In this article, search regions are defined by at most two kinematic cuts 
applied after predefined selection criteria. These selection criteria, motivated by typical 
LHC searches for jets and missing energy in all-hadronic channels, are: 

• a jet pt selection of pxi > 100 GeV and p-ri > 50 GeV, for 2 < i < channel multiplicity; 

• a pti > 20 GeV veto on all isolated leptons; 

• a \ A<f)(ji, $t)\ > 0.2 cut between that transverse momentum vector of the i th jet and 
the missing transverse momentum vector, where i = 1,2,3 2 ; 

• a \r]j t \ < 2.5 requirement on the momentum vector of the i th jet, where i — 1,2, 3 3 ; 

• a $t/M s s > 0.30, 0.25, 0.20 cut for dijet, trijet and multijet channels, respectively, for 
controlling QCD backgrounds. 

The expected 2<r-sensitivity quantifies how well a search region performs for a simplified 
model at a given point in the nig, m x o and cr pro d parameter space. A given model is declared 
to be within the 2<r-reach of a search region if the number of signal events passing the search 
region cuts is equal or greater than twice the background uncertainties, both statistical and 
systematic: 



iV S i gn ai > A£? 2CT = 2 x ^(AB stat y + (A5 syst ) 2 , (16) 

where a 30% systematic background uncertainty, A£> syst = 0.30iVbk g , is assumed throughout 
this article. The statistical error on the number of background events, Ai? sta t, is the Poisson 
fluctuation of the expected number of background events and is equal to >jN hkg in the 
Gaussian limit [63J. For a given search region, the minimum value of a(pp — > gg) x B(g — > 



nj + $ T ) that satisfies relation (16) can be computed 

AZ? 2<J 

a% od (pp ggX) x B(~g nj + $ T ) = AxtxC (17) 

where A is the acceptance and e is the efficiency of the search region's cuts for the signal. 

The reach of the LHC depends on the choice of search regions. In order to extract 
the optimal reach for each point in parameter space, a large ensemble of search regions in 



1 M e s is defined as M c g = $t + ^2 Pru where n is the channel's jet multiplicity. 

i=l 

2 For the dijet channel, this cut is only applied to the leading two jets. 

3 For the dijet and trijet channels, this cut is only applied to the leading two or three jets, respectively. 
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the kinematic variables discussed above was explored. In the final outcome, it turned out 
that combined cuts in {$t-> Ht}, optimized individually for each simplified model, yielded 
the best sensitivity to the production cross section times branching ratio of the parameter 
space explored. In practice, that was obtained by simulating simplified models in a grid of 
{rrig, m x o} for all six decay modes considered, in 50 GeV steps for m g and 25 GeV steps for 
m x o. Each simulated simplified model had its reach evaluated in a grid of combined cuts in 
{$ T ,H T }, with 100 GeV < $ T < 600 GeV in 50 GeV steps and $ T < H T < 1200 GeV in 
50 GeV steps. 

The optimal sensitivity on the cross section times branching ratio for all mass spectra 
and decay modes is displayed in Fig. [3] and App. [A} 



1-step cascade decay (r=1/2) 1-step cascade decay (r=1/2) 




m g (GeV) m g (GeV) 



FIG. 3: Contours of the 2<7-sensitivity for a(pp gg) x B(g —?■ qq'x^ —> qq'iW^x)), with 
C = lib" 1 at 7 TeV. The expression corresponds to the g pair-production cross section times the 
branching ratio for the one step cascade decay mode, where the intermediate particle mass lies at 
m x ± = m x o + l/2(m g — m x o). The optimal reach (right) is contrasted with the reach provided 
by the 1.04 fb _1 ATLAS search for supersymmetry in jets plus missing energy presented at EPS 
2011 [57] (left). The contour values are specified on the color scale on the right. The red lines 
delimit the expected reach of models for which <7 pro d x B g is a simple parametrization of the g 
NLO-QCD production. They correspond to: (i) cr pro d x Bj = 3 x 0NLO-QCD (dashed line), (ii) 
<Vod x &g = ^nlo-QCD (solid line), and (iii) a prod x B g = 0.3 x onlo-QCD (dotted line). The blue 
lines in the left plot are analogous to the red lines, but correspond to the estimated 2cr-exclusion 
regions from the data. 

Fig. [3] contrasts the optimal 2cr-sensitivities at 7 TeV with the reach of the 1.04 fb _1 
ATLAS search for supersymmetry in jets plus missing energy presented at EPS 2011 |57J. 
The contours specify the minimum accessible values of the g production cross section times 
branching ratio, cr pro d x B g , as a function of nig and m x o, for the one-step cascade decay 
with r = 1/2. Analogous plots displaying the optimal reach for the other decay topologies 
considered in this paper are contained in App. [Ajfor two luminosity scenarios, C = 40pb _1 
andlfV 1 . 
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The search regions used in the analysis in [57] are schematically displayed the table below: 



Search Region of [57] 


> 2 jets 


> 3 jets 


> 4 jets 


High Mass 


$ T { GeV) 


> 130 


> 130 


> 130 


> 130 


PTi(GeV) 


> 130 


> 130 


> 130 


> 130 


p T 2( GeV) 


> 40 


> 40 


> 40 


> 80 


p T3 (GeV) 


— 


> 40 


> 40 


> 80 


PT 4 (GeV) 






> 40 


> 80 


A0(p T , $ T ) 


> 0.4 


> 0.4 


> 0.4 


> 0.4 


$t/M cS 


> 0.3 


> 0.25 


> 0.25 


> 0.2 


M eff (GeV) 


> 1000 


> 1000 


> 500/1000 


> 1100 


Background 


62 3 ±4 3 

UZ.O ±g 2 


c,e;±3.8 
od ±7.3 


984±3Sb/33.4±« 


13 2 ±L9 

±0 - z ±2.6 


Data 


58 


59 


1118/40 


18 


excluded e x cr(fb) 


24 


30 


477/32 
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Fig. [3] also shows the difference in mass reach when optimized cuts are used versus the 
cuts used in [57] for a few rescalings of cr pro d x B g in terms of the g NLO-QCD production 
cross section. Given that for the analysis of [57] the actual data and background estimates 
were available, both the expected 2<r-sensitivity and the estimated 2a-exclusion regions are 
displayed (red and blue lines, respectively). The difference between expected and estimated 
exclusion regions is due to a 1-sigma excess in the measured number of events in the multijet 
signal regions relative to estimated backgrounds. 

Inspection of the sensitivity contours in Fig. [3] reveals that, although the ATLAS search 
regions [57j were well optimized for heavy spectra and large mass splittings, it compromised 
its sensitivity for spectra with nig < 700 GeV and m x o < 300 GeV. One may naively argue 
that this region has already been excluded and therefore there is no reason for optimizing 
search sensitivities to nig < 700 GeV. That statement is true, however, under the hypothesis 
that the gluino production cross section is purely determined by its QCD interactions and 
that its branching ratio for this 1-step cascade decay mode is equal to 100%. As discussed 
in Sec. [2j this may not be the case. The production cross section for the gluino may be 
suppressed relative to the reference QCD cross section through t-channel squark interference, 
or the gluino may have a significant branching ratio to invisible or soft modes, in which case 
these lighter spectra remain experimentally viable. 

Given that there is no particular reason to expect that CMSSM type-spectra are more 
likely than lighter or more compressed spectra, it is clear that the LHC discovery potential 
on jets plus MET searches may be crippled by a limited choice of search regions. On the 
other hand, measurements of Standard Model backgrounds is a challenging task that may 
render impractical those analyses demanding several dozens of search regions. The viability 
of a search strategy that is sensitive to all space of signatures will be addressed in the next 
section, where a thorough investigation of the $t and Ht cut-space will be performed, and 
guidelines for a comprehensive search strategy will be proposed. 
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5. COMPREHENSIVE SEARCH STRATEGIES 

This section addresses the challenge of creating effective and comprehensive search strate- 
gies, which is how to discover any theory from the entire space of models with the minimum 
amount of integrated luminosity and a practical number of search regions. There is no single 
search region whose coverage is close to optimal in the whole parameter space of theories. 
On the other hand, there is a practical limitation to the number of search regions that can 
be studied on a given analysis, due mainly to the challenge posed by measuring the Standard 
Model background in each search region. A useful approach to this problem that has been 
regularly implemented is to design searches that optimize their sensitivity to the regions of 
parameter space that are theoretically motivated. Of course, theoretical motivation is sub- 
jective and time dependent. This study adopts a different perspective, namely, to remain 
agnostic about the likelihood of a given range of model parameters and look for a minimal, 
comprehensive set of search regions whose combined reach is close to optimal in the full 
model space. This approach will be referred to as a Multiregion Search Strategy in jets plus 
missing energy at the LHC. 




trig (GeV) m g (GeV) m t (GeV) 



FIG. 4: Reach of the search region S = {$ T > 150 GeV, H T > 750 GeV} in the inclusive multijet 
channel for the six different decay topologies. The efficacy of this search region is £ < 1.10 in the 
dark purple region, £ < 1.20 in the medium purple region and £ < 1.30 in the light purple region, 
assuming an integrated luminosity of lfb -1 . 

As mentioned in the previous section, an extensive set of search regions were considered 
in order quantify the optimal sensitivity to our simplified models, characterized by cuts on 
missing and visible energy within the following range : 100 GeV < $ T < 600 GeV and $ T < 
Ht < 1200 GeV. Alternative two-cut search regions, e.g. $t and M e g, were investigated as 
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well but were not found to be superior to $t and Ht- Each search region is efficient in a 
portion of the model parameter space. Conversely, each point in the parameter space may 
be covered by several search regions. The goal is to find a minimal set of search regions that 
covers the entire model parameter space. 

Efficacy, £, is a variable that quantifies how close to optimal the reach of a given search 
region S is. Efficacy is defined as 



a prod (S;M) 



where M. denotes a specific simplified model with a given mass spectrum and decay topology, 
S is the given search region (specified by a set of cuts), cr pro( i(<S; M) is the estimated reach 
in the production cross section for M. using the search S. cr°^ d (A^) is the optimal estimated 
reach in the production cross section for M. when considering all search regions, 

<r*UM) = min {a prod (S; M)\S} . (19) 

By definition, £ is greater than or equal to unity for any search region. The closer S(S, M.) is 
to unity, the more sensitive S is to M.. Choosing search regions that have E ~ 1 decreases the 
amount of integrated luminosity necessary to discover new physics. The goal of this article 
is to find a set of search regions that covers all possible models with £ < £ C rit- To accomplish 
this, £ (S; M.) must be computed over the entire space of S and Ai. As an illustration, Fig. [4] 
shows the reach of the search region S = {$t > 150 GeV, Ht > 750 GeV} in the inclusive 
tetrajet channel for efficacies £ crit = 1.1, £ CT it = 1.2 and £ crit = 1.3. The next step is to find 
a minimal set of search regions whose combined reach spans all the simplified models and 
decay topologies under consideration. 



5.1. Multiregion Search Strategy in $t & Ht 

The search region optimization problem is computationally intensive and has no unique 
solution. The approach adopted in this study was to perform the optimization with a 
genetic algorithm. A random "population" of multiple search regions is initially generated. 
Mutations and performance selection of the elements of the population are successively 
implemented, until one or more multiple search regions with the desired criteria are found. 

One of the multisearch regions found by our genetic algorithm with a global efficacy 
requirement of £ cr n = 1.3 is displayed in Fig. |5j The optimization was performed for three 
simultaneous values of integrated luminosity: C = 10 pb -1 , 100 pb -1 , 1 fb~ . It was found 
that six search regions were necessary to cover the full model space with £ cr ; t = 1.3. The 
search regions for this particular solution are: 
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where the background cross section is broken down into its main contributions ( W +j, Z° + 
nj, ti+j). The QCD contribution to drops rapidly. At > 100 GeV, the cross section 
is 0(120 fb) and by $t > 150 GeV the cross section has dropped to 0(3.5 fb). Only search 
region 3 had any appreciable QCD contribution to the background and was estimated to be 
~ 3.7fb. Detector effects contributing to px, modeled by PGS4, have large uncertainties, 
frequently arising from reducible backgrounds. Therefore the QCD contribution to these 
search regions is not quoted, but is believed to be subdominant to the listed backgrounds 
in all search regions. 

The minimum number of search regions depends on the global optimality requirement, 
i.e., on £ cv \t- For £ crit = 1.2, the number of search regions increases from 6 to 13. If £ cvlt is 
relaxed to 1.5, a set of 4 search regions is sufficient. Each search region captures a specific 
kinematic regime, although the values of the $t and cuts have some room for variation 
without affecting the coverage. The qualitative kinematic regimes of each search region and 
the patches of phase space covered by each one of them are described below. 

• Dijet high MET 

A high $t cut on the inclusive dijet channel is required to provide coverage of theories 
with nearly degenerate spectra [22J, [23j El] . This qualitative region is expected because 
compressed spectra are only visible in g events that come with initial or final state 
radiation. The dijet channel was also important in the heavy g - light x° region for 
the two-body direct decay mode, where the events typically have two hard jets and 
very energetic x°' s - 

• Trijet high MET 

The trijet inclusive channel is important for providing sensitivity to heavy g's directly 
decaying via 2-body as well. In contrast to the dijet channel, the trijet is more sensitive 
to slightly heavier x° masses. 

• Multijet low MET 

All the other regions were best covered by cuts on the inclusive tetrajet channel. A 
search region with a low $t cut is necessary to provide sensitivity to light g's, specially 
for the decays that yield high jet multiplicities. 

• Multijet high H T 

A high visible energy cut was essential for coverage of models with heavy g's and large 
mass splittings in their spectra. For those types of models, the heavy g decays yield 
events where a large fraction of the energy going into visible states, allowing the cuts 
to go far in the distribution tails where backgrounds die off rapidly. 

• Multijet moderate MET 

A moderate $t cut is useful for intermediate trig in regions of moderate mass splittings. 
It is also important for light g's decaying via 2-body because this decay mode typically 
yields fewer but harder jets and more missing energy in comparison to 3-body or 
cascade decays. 

• Multijet high MET 
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FIG. 5: Minimal multiple search region on $t and Ht whose combined reach is within 30% of 
optimal for all kinematical regions and decay topologies for the integrated luminosity range C = 
10pb -1 — lfb -1 . The dark dots correspond to benchmark simplified models that are representative 



of the full phase-space and can be used in optimizing searches, see Tab. 3. A 3.F 



Finally, a high search region plays an important role for models with heavy ^'s 
decaying to non-relativistic x°' s - I n these types of spectra, x° carries off a large fraction 
of the energy available in the decay, resulting in lower H T but still a sizable amount 
of $ T . 

Fig. [6] illustrates the complementarity of the multiple search regions and the limitations 
of an isolated search region. It displays the Ht distribution in the multijet channel for all 
the dominant backgrounds and two signal points: 

• trig = 200 GeV, with a 25% branching ratio into the 1-step cascade decay mode with 
r = 1/2 and m x o = 55 GeV and the remaining decays invisible, 

• m g = 800 GeV, with a 100% branching ratio into the 1-step cascade decay mode with 
r = 1/2 and m x o = 55 GeV. 
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FIG. 6: Ht distribution in the multijet channel with a $t > 100 GeV requirement. The dominant 
backgrounds, tt + jets and W/Z + jets, are displayed in gray, whereas to signal points are displayed 
in color. The Multijet Low MET search region, with a Ht ^ 450 GeV cut, would be sensitive to 
the 200 GeV g signal, whereas the Multijet High Ht search region would be sensitive to the 800 
GeV g signal. 



The Multijet Low MET search region would cut on Ht around ~ 450 GeV, making the 
rrig = 200 GeV signal discoverable, but swamping the m g = 800 GeV signal in backgrounds. 
In order to make the latter signal visible, a hard H T cut, > 1000 GeV, is required, as 
provided by a Multijet High Ht search region, which would however completely kill the 
nig = 200 GeV signal. 

This Multiregion Search Strategy is crucial to ensure that new physics signatures will 
not be missed because there have been relatively few studies of the full signature space. 
Performing a measurement in a single search region limits the reach to a preferred corner 
of model space. Only the combination of complementary search regions captures the full 
potential of the LHC, providing near-optimal coverage. 

Fully studying any simplified model requires extensive Monte Carlo calculations of the 
different signals. The primary concern of this article is to ensure that no signal is missed 
and to motivate more thorough searches in jets and missing energy. Simulation of hypothet- 
ical signals is computationally costly and minimizing Monte Carlo generation is desirable, 
particularly for the experiments at the LHC where full detector simulations limit the size of 
surveys of hypothetical signals. The search regions found here need to be reanalyzed by the 
experimentalists to make sure that they do provide the coverage claimed and are not limited 
by unforeseen backgrounds. In order to facilitate this process, a selection of a few dozen 
benchmark simplified models is provided to ensure that sensitivity is not lost when tuning 
searches. These benchmark simplified models are fully specified in App. [Bj Tabs. 3. A 3.F 
and marked as dark dots in Fig. [5j 
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5.2. Alternative p? Selection Criteria 

The multiregion search strategy in $t and presented in the last subsection used 
the pre-selection criteria described in Sec. |4j These pre-selection criteria are determined 
primarily by triggering and can in principle be tightened to provide better limits. For the 
three-body direct decay, two-step cascade decay and r = |, | one-step cascade decay modes, 
30% to 40% sensitivity can be gained in the range 350 GeV < rrig < 600 GeV, m x o < 150 GeV 
by hardening the p^-selection criteria on the sub-leading jets of the inclusive tetrajet channel: 

Pti, Pt2, Pts, Pt4 > 100 GeV. (20) 

These selection criteria could replace the ones in Sec. |4j with a gain in sensitivity in the 
low m x o region for some decay topologies. However, this tighter pre-selection results in a 
loss in sensitivity ranging between 20% to 40% in the whole parameter space for the r = § 
one-step cascade decay mode. It also renders the tetrajet channel inefficient for the two- 
body direct decay mode. Tightening the pre-selection criteria in lower multiplicity channels 
to pt > 100 GeV results in a 10% to 40% loss of sensitivity for more compressed spectra 
and would universally affect the reach for m g — m x o < 300 GeV. This discussion shows that 
the search regions defined in the previous subsection can be improved, but at the cost of a 
significantly more complicated design. 

6. MULTIPLE DECAY MODES 

The studies considered so far are applicable only when a single g decay mode contributes 
to the all-hadronic jets plus missing energy channel. When the g has two or more decay 
modes contributing to this channel, translating results from the single g decay mode into 
these more generic decay patterns is not completely straightforward. This section addresses 
how to infer sensitivities of searches to simplified models when there are multiple decay 
modes, but only single decay modes have been explicitly studied. Understanding how to use 
models with single decay modes in a more general context greatly enhances the applicability 
of simplified models and reduces the overall complexity of studying more general examples. 

In the following we provide a quantitative illustration of an estimation of the sensitivity 
to the cross section for models with multiple g decay modes. Consider a model in which g 
decays through two modes, A and B, with branching ratios Ba and respectively. The 
goal is to be able to set a limit on the g pair-production cross section with these two decay 
modes, 

Oprod(#4, B B \S), 

when only the quantities 

Oprod(l,0|S), CTprod(0, 1|<S) 

are known. 

A conservative limit on the production cross-section can be extracted by considering the 
most constrained decay mode and ignoring the other decay modes, 

in ■ f°"prod(l,0) CTprod(0,l)\ , 

a prod {B Al B B ) < mm j — 2 , — 2 j. (21) 
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This can be a considerable underestimate of the actual search sensitivities because it assumes 
that the mixed and less sensitive decay modes make no contribution to the signal region. 

The sensitivity of a given search region S depends on the number of events of each 
topology that pass the cuts defining S. The topologies AA, BB and AB contribute the 
following number of events to the signal region 



N A a(S) = e AA (S)B 2 A x N 
N B b(S) = t BB {S)B 2 B x N 
N AB (S) = 2e AB (S)B A B B x N 



for gg — > AA, 
for gg — > BB, 
for gg — > AB, 



(22) 



where e(S) is the efficiency for the given decay mode to pass the cuts defining the search 
region S and 

N(S) = Cxa prod (B A ,B B \S) 
is the total number of events. It is useful to define the weighted efficiency for S 

e T (B A , B B \S) = e AA (S)B\ + e BB (S)J3 2 B + 2t AB {S)B A B B . (23) 
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FIG. 7: Benchmark model nig = 700 GeV, m x o = 80 GeV with two g decay modes: A = two- 
body direct decay and B = two-step cascade decay with r = r' = |. Three topologies are possible: 
non-hybrid AA, BB and hybrid AB. This example is studied in the text with B A = B B = 50%. 



As a specific example, Fig. [7] illustrates a theory with benchmark masses 



700 GeV, 



m x° 



80 GeV in the case in which mode A corresponds to the 2-body di- 



rect decay of g and mode B to the 2-step cascade decay. Tab. [T] displays the efficiencies and 



background estimates for the search regions of Sec. 5 5.1 



The expected 2<r-sensitivity to the cross section is extracted by demanding that the 
number of signal events in S, N(S) = Gt(S) x N, be greater than twice the background 



uncertainty. From (17), the 2a limit on the production cross section from S is 



a£ od (B A ,B B \S) 



AB 2a (S) 



Cxe T (B A ,B B ,S) 



(24) 
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where AB 2a is the 2a uncertainty in the background given in (16). 

The question of how to extract limits for theories with multiple decay modes reduces to 
how much is known about e AB (S) when only e AA (S) and e BB (S) are known. Although it is 
possible to obtain the efficiencies for hybrid events through separate Monte Carlo studies, 
when there are several different decay modes it is computationally expensive to extract 
the efficiencies for all hybrid events (\Nd(Nd + 1) additional Monte Carlo studies when Nj, 
different decays are present). A more conservative approach is to estimate the efficiencies 
for hybrid events using the known efficiencies for non-hybrid events. If the decay modes A 
and B are relatively similar, then the efficiencies for the hybrid events are bounded by 

min {e AA (S), e BB (S)} < e AB (S) < max {e AA (S), e BB (S)} . (25) 

When this two-sided bound is satisfied, a pio ^(ggX; B A , B B \S) can be bounded from both 
sides 



<&d(M|S) 2CT m n lq)< < od (i,o|s) 

B\ + (B 2 B + 2B B B A )^ ~ a P rodl ° A ' ° B]0) ~ B\ + 2B A B B + B 2 B ^ 



(26) 



assuming e AA < e BB . 

This two-sided bound on efficiencies for hybrid decays is not always satisfied when the 
decay modes are substantially different, e.g. two-body direct decay vs. two-step cascade 
decay. In general there is no way to bound e AB in terms of e AA and e BB . However, for the 
search regions studied in this article, the efficiencies for the hybrid decays do satisfy a bound 
of 

e AB (S) > mm{e AA (S),e BB (S)} , (27) 

which was verified empirically. By taking the lower bound of the relation above, a conser- 
vative estimation of the cross section sensitivity can be placed, 



B\ + 2B B B A ^- A + B|g* - ^ d ^' bBlS) ~B\ + 2B A B B + B^ ^ ] 

assuming that e AA < e BB . The lower limit on a 2 ° od {B Al B B \S) arises from the constraint 
€ AB < 1 and is primarily useful when the AB topology branching ratio is subdominant to 
the AA or BB branching ratio. In Tab. [2] this conservative estimate of the cross section is 
applied to the specific example of Fig. [7j with B A = B B = 50%. Notice that even though the 
search regions were not explicitly designed to be sensitive to these hybrid decays, sensitivity 
is nearly optimal with S = 1.02 for the most sensitive search region. The optimal sensitivity 
for this benchmark simplified model is set by the high-if^ cut of the multiple search region. 
The conservative estimate for the limit on the production cross-section (162 fb) is higher 
than the actual value (108 fb) by 49%, and higher than the optimal value (106 fb) by 
~ 53%. More generally, in the framework of comprehensive multiregion search strategies, 
this conservative estimate on the cross section limit is within 50% of the actual limit for all 
of the simplified models studied. 



7. DISCUSSION 

This work focused on the optimization of searches for new colored states with jets plus 
missing energy signatures. Simplified Models were used to capture relevant new physics 
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450 
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0.25 


0.20 


0.43 


4 


4+i 


150 


950 


13.4 


10.4 


0.075 


0.064 


0.13 


5 


4+i 


250 


300 


105.7 


66.6 


0.21 


0.078 


0.33 


6 


4+J 


350 


600 


18.1 


13.8 


0.14 


0.017 


0.15 



TABLE 1: Signal efficiencies for the multiple search region of Sec. 5 5.1 for benchmark masses 
rrig = 700GeV, m x o = 80GeV and g decay modes A = 2-body direct decay and B = 2-step cascade 
decay. Also included are the expected background cross section in the signal region <7bkg X e and 
corresponding statistical and systematic uncertainties AB 2a for a luminosity of 1 fb _1 . 



s 


Cd(M) 


^prod(0,l) 


^ctual(0-5,0.5) 


a c 2 o CT ns .(0.5,0.5) 


rr 2a ln 2a 
"cons. / "actual 


^actual 


£cons. 


1 


101 fb 


6690 fb 


267 fb 


389 fb 


1.46 


2.52 


3.67 


2 


95.3 fb 


4260 fb 


216 fb 


357 fb 


1.65 


2.04 


3.37 


3 


1610 fb 


2020 fb 


1234 fb 


1901 fb 


1.54 


11.6 


17.9 


4 


144 fb 


169 fb 


108 fb 


162 fb 


1.49 


1.02 


1.53 


5 


317 fb 


854 fb 


281 fb 


600 fb 


2.13 


2.65 


5.66 


6 


98.6 fb 


812 fb 


121 fb 


289 fb 


2.39 


1.14 


2.73 



TABLE 2: Cross section sensitivity for the benchmark masses rrig = 700GeV, m x o = 80GeV and 
the two g decay modes A = 2-body direct decay and B = 2-step cascade decay, where, in order to 
maximize the number of hybrid events, we take Ba = = 50%. The highlighted orange search is 
most effective for the A A topology and the yellow is most sensitive for the BB and AB topologies. 
The actual sensitivity (T ac t U ai was computed using (24) and the efficiencies displayed in Tab.Jl} The 
conservative estimate a cons in case the efficiency for hybrid events is unknown is obtained by taking 



the lower bound in (27) for €ab- The last column displays the efficacy of each search region under 



the conservative estimates ct coi1s . The efficacy, £, is defined in (18) and quantifies how close the 
cross section limits are from the optimal one, <7 op tmiai = 106 fb. 



features with the simplest spectra in the most important decay topologies. With optimized 
search regions, the reach of the 7 TeV LHC for heavy colored octets decaying to jets and 
a long-lived invisible particle was estimated in a wide range of masses and the following 
decay topologies: two- and three-body direct decays, one-step cascade decays and two-step 
cascade decays. 

The optimal reach requires tuning the search regions for each individual spectrum, and 
is therefore not practical. A more minimal search strategy, nearly as effective as a fully 
optimized search, was presented, consisting of a set of six search regions. These search 
regions are characterized by cuts on missing and visible energy, whose combined reach is 
within 30% of optimal for all kinematic regimes, decay topologies, and integrated luminosities 
in the range 10 pb" 1 < £ < lfb -1 . Although not unique, multiregion search strategies 
share qualitative features that capture specific regions of the phase space of signatures. For 
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instance, a hard missing transverse energy cut in the inclusive dijet channel is required for 
coverage of compressed spectra. Other regions are best covered by higher jet multiplicities, 
in particular the tetrajet channel. 

The efficacy of the search regions depends on the assumptions about the backgrounds, in 
particular the systematic uncertainty. For instance, the studies in this article used a 30% sys- 
tematic uncertainty on all backgrounds. If there are larger systematic uncertainties, search 
regions with harder cuts will be preferred so as to remove backgrounds, whereas looser cuts 
will be preferred for smaller systematics. Since these uncertainties are moving targets, a set 
of benchmark simplified models are given in App. |B]that provide a representative sampling 
of the whole space of models and can be used for designing a comprehensive multiregion 
search strategy. These benchmark models can be used to optimize search strategies with 
more realistic background calculations or background measurements. 

The results of this work reveal a promising picture for the upcoming LHC results that 
will shortly be released. Specifically, with 45 pb -1 , the LHC will be able to test g masses 
up to 600 GeV for light x°' s > an d have nearly complete coverage, independent of m x o, up to 
350 GeV. With lfb" 1 , the reach on g masses will extend up to 850 GeV and have complete 
coverage up to 400 GeV. Jets and missing energy are the first channels to look for new 
physics with QCD interactions below the TeV scale, and multiregion search strategies will 
play a key role in the discovery process. 
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Appendix A: Reach Estimates 




FIG. 8: Contours of the 2cr-sensitivity for a(pp — > gg) x B for the six decay modes considered in 
this work, assuming C = 45 pb -1 and *Js = 7TeV. The contour values are specified on the color 
scale on the left. The red (dashed) lines are simple parametrizations of <7p ro d x B in terms of the g 
NLO-QCD production corresponding to <7 pro d x 6/<tnlo-QCD = 3.0, 1.0,0.3, moving from right to 
left, respectively. 




FIG. 9: Same as in Fig. M but with 1 fb 1 of integrated luminosity at y/s = 7 TeV. 
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Appendix B: Benchmark Models 

This appendix lists the fully specified benchmark models. 



Name 


m- g (GeV) 


m x o (GeV) 


^fpb-i (P b ) 


(Pb) 


Cd (P b ) 


gr 


350 


150 


13 


2.9 


24 




550 


300 


3.2 


0.69 


1.2 




800 


100 


0.23 


0.049 


0.061 




800 


350 


0.46 


0.097 


0.061 


gr 


850 


840 


14 


3.3 


0.036 


gr 


900 


190 


0.17 


0.036 


0.021 



TABLE 3.A: Benchmark simplified models for <?'s that decay directly to x through a two-body 
decay. The optimal reach in the cross section for each one of the benchmark models, a opt , is 
displayed for two luminosity scenarios: 45 pb -1 and lfb -1 . Also displayed is the reference NLO- 
QCD cross section for g pair-production. 



Name 


m- g (GeV) 


rriyX) (GeV) 


^spb- 1 (P b ) 


(Pb) 


<od (Pb) 




250 


50 


32 


9.5 


180 




400 


100 


5.7 


1.2 


14 




400 


350 


23 


5.9 


14 


gr 


650 


300 


1.2 


0.26 


0.34 



TABLE 3.B: Benchmark simplified models for gr's that decay directly to x through a three-body 
decay. 



Name 


m g (GeV) 


m x ± (GeV) 


m x o (GeV) 


^spb- 1 (Pb) 


(Pb) 


< C od (Pb) 


gl cv 


150 


75 


50 


214 


156 


2900 


gr 


400 


160 


80 


10 


2.3 


14 


gr 


450 


375 


350 


17 


4.2 


4.8 


gr 


600 


300 


200 


2.2 


0.48 


0.62 


gr 


850 


272.5 


80 


0.30 


0.064 


0.036 



TABLE 3.C: Benchmark simplified models for g's that decay through a one-step cascade to 
with r = 1/4, where r = (m x ± — m x o)/(m,g — m x o). 
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Name 


m 5 (GeV) 


m^± (GeV) 


m^o (GeV) 


f?? 1 , 1 (pb) 

45 pb -1 VK ; 


c?^* i (pb) 

lib -1 1 


a QC » (pb) 

prod V- 1 ^ / 


r;iCD 




ZZO 


inn 
ZUU 


57 


zo 


loU 


or 


300 


175 


50 


27 


7.4 


62 


C 1CD 


550 


525 


500 


18 


4.5 


1.2 


C 1CD 
»9 


700 


550 


200 


0.84 


0.18 


0.19 


C 1CD 
» 10 


900 


500 


100 


0.20 


0.042 


0.021 



TABLE 3.D: Benchmark simplified models for g's that decay through a one-step cascade to %' 
with r = 1/2. 



Name 


m § (GeV) 


m x ± (GeV) 


m x o (GeV) 


^spb-i (P b ) 


^ift-i (Pb) 


<od (Pb) 


C 1CD 

»n 


250 


187.5 





37 


12 


180 


C 1CD 
» 12 


350 


312.5 


200 


27 


7.5 


24 


C 1CD 
»13 


450 


362.5 


100 


4.4 


0.94 


4.8 


C 1CD 
»14 


850 


837.5 


800 


15 


3.4 


0.036 


C 1CD 
f 15 


900 


775 


400 


0.42 


0.090 


0.021 



TABLE 3.E: Benchmark simplified models for g's that decay through a one-step cascade to %' 
with r = 3/4. 



Name 


m § (GeV) 


m x ± (GeV) 


m x /o (GeV) 


m x o (GeV) 


<5 P b-i (Pb) 


(pb) 


Cd (Pb) 


or 


300 


175 


112.5 


50 


37 


12 


62 


Qr 


750 


450 


300 


150 


0.67 


0.14 


0.11 


or 


750 


650 


600 


550 


6.6 


1.4 


0.11 


or 


800 


775 


762.5 


750 


15 


3.4 


0.061 


or 


900 


450 


225 





0.33 


0.070 


0.021 



TABLE 3.F: Benchmark simplified models for g's that decay through a 2-step cascade to x with 
r = r' = 1/2, where r = (m x ± — m x o)/(m g — m x o) and r' = (m x io — m x o) / (m x ± — m x o). 
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